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The non-perturbative quark-gluon interaction depends significantly on the value of the quark 
mass. In particular, in the light quark mass regime, correlations are strongly influenced by dy- 
namics associated to chiral symmetry breaking. We use the Interacting Instanton Liquid Model 
(IILM) as a tool to investigate the microscopic dynamical mechanisms which underly the depen- 
dence on the quark mass and drive the transition into the chiral regime of QCD. To ensure the 
validity of the model, we first verify that the dependence on the quark mass for several observ- 
ables calculated in the IILM agrees well with the predictions of chiral perturbation theory and with 
lattice simulations. We then show that a quark mass m* sw 80 MeV emerging naturally from the 
model specifies the mass scale above which the dynamics associated with low-lying eigenmodes 
of the Direac operator becomes sub-leading and the contribution of the fermion determinant is 
suppressed. 
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1. Introduction 

Although Lattice QCD provides the only available ab-initio tool for solving QCD non-perturbatively, 
such a framework does not allow a direct identification of the dynamical mechanism underlying 
hadron structure and of the relevant degrees of freedom at different scales. In this talk we present 
our recent attempt to extract complementary insight using the Interacting Instanton Liquid Model 
(IILM) [1]. 

In the IILM, the QCD path integral over all possible gluon field configurations is replaced by 
an effective theory in which instantons are the effective degrees of freedom, and the gauge fields 
of the theory are those generated by integrating over the positions, color orientations, and sizes of 
instantons. The QCD partition function reduced to: 

i f N++N- N f 

~¥qcd — 2?ilm = 52 ~xj ur r / ' dn i J(p i )e- Sw r[ det (^ + ™/)- C 1 - 1 ) 

N + ,N- N +- N - - J i i 

Here, d£2, = dt/,-d 4 z,dpj is the measure in the space of collective coordinates, color orientation, posi- 
tion and size, associated with the single instantons. Quantum fluctuations are included in Gaussian 
approximation, through the semi-classical instanton amplitude d{pi). Si nt is a bosonic interaction 
between pseudo-particles which includes a phenomenological short-range repulsive core required 
to remove large-sized instantons from the vacuum. In the formulation of the model we have con- 
sidered, the strength of such a repulsion is the only phenomenological parameter, which has to be 
tuned to reproduce observations. 

Numerical simulations in the Instanton Liquid Model performed in the last few years have been 
shown to be in very good agreement with both phenomenology [2, 3] and lattice simulations [4]. 
Here, we would like to understand the mass range in which instanton mediated chiral dynamics 
is manifest and the extent to which it is described by chiral perturbation theory. In particular, 
we would like to identify the mass scale or scales above which the contribution of the fermion 
determinantis suppressed, and the near-zero modes of the Dirac operator become sub-dominant. 

2. The Chiral regime in the IILM 

In order to investigate hadron structure in the chiral regime of QCD we first check that IILM 
gives predictions in agreement with #pt. To test the IILM in the chiral regime, we have analyzed 
the structure of the spectrum of the Dirac operator. In ^pt there are clear statements on how the 
distribution of eigenvalues of the Dirac operator p(A) changes with the mass of current quarks. In 
the chiral limit, m q = 0, the Banks-Casher relation implies that p(A) becomes flat at the origin [5] 

lim lim p(A) =Const+0(A 2 ) (N f = 2), (2.1) 

In Fig. 1 we present the Dirac spectrum computed in the IILM for m q = and Nf = 2. Clearly, in 
the X —y limit, we observe that the spectrum develops the expected plateau. 

If the quark mass is small but finite, #pt predicts a peak near the origin of the form p(A) — 
p(0) oc m a \ x ~ a . Such a particular structure of the Dirac spectrum away from the chiral limit is 
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A [Aqcd] 

Figure 1: Spectral density of the Dirac operator for N c = 3, Nf = 2 and m q = 0. 



reflected in a specific dependence of the quark condensate on the quark mass: 

(qq)=-2m q TdA-^z ( 2 - 2 ) 
Jo m l q + A, z 

Using this relationship, we can obtain the IILM prediction for the condensate for different m q , 
which can be compared with the &{p A ) #pt formula 

We have found that the observed quark mass dependence is consistent with the form expected from 
chiral perturbation theory. In addition, the low-energy constants /o, Bo and Z3 extracted from a 
best-fit of the IILM data using (2.3) are consistent with the value extracted from phenomenology 
and recent lattice results [6] : 

fo = 0.085 ±0.005, S = 2.4 ±0.2, 7 3 = 3.80 ±0.04, (2.4) 

which give (0|^|0) = (0.26 ±0.01 GeV) 3 . 

3. Pion and Nucleon mass in the IILM 

The second step in our investigation of the instanton contribution to the low-energy sector 
of QCD is to verify that the IILM can be used to calculate hadron masses and that these masses 
depend on the quark mass in the specific way predicted from #pt. We have considered pion and 
nucleon because these are the two hadrons in which the instanton-induced interaction is most in- 
tense. To calculate the hadron masses we have adopted standard effective mass plot analysis, which 
is routinely used to extract lowest-lying hadron masses from lattice simulations. 

To compute the mass of pion and nucleon we have evaluated the correlation function: 

G*(t) = J d 3 x(0|r[7 5 feT)7 5 (0,0)|0), G N (z) = I d 3 x(0\T[j N (x, z)j N (0,0)P+\0), 
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Figure 2: Dependence of the quark condensate from the quark mass, the fit is obtained using the (2.3). 



where the two currents are given by j%{x) = £ u a (x)CYsd {x)u c (x) and jf(x) = q{x)y^z a q{x), 
and P + = is the positive-parity projector. 

The mass of these hadrons can then be extracted from the plateau in the large Euclidean time 
limit of the effective mass, i.e. using 



We have computed these correlation functions for five different values of the quark mass in 
the range 20 < m q < 90 MeV, finding the results reported in Table 1 . We have found that all the 
pion effective mass plots display a very clean plateau, from which it is possible to unambiguously 
read-off the pion mass, from a correlated chi-square fit. In the interpolation of the effective mass 
plots of the nucleon not only statistical but also systematic error has to be taken into account. 

Table 1: Pion and nucleon masses fitted by effective mass plot with % 2 jndf < 1. The quark masses are 
determined at a scale 2 GeV as in [7, 8] 



m q [GeV] Pion [GeV] Nucleon [GeV] 



0.021 
0.03 
0.05 
0.07 
0.09 



0.300 ±0.004 
0.360 ±0.004 
0.460 ±0.004 
0.530 ±0.004 
0.600 ±0.004 



+0.05 
-0.05 
+0.01 
-0.07 

1 20+ 002 
1 - zu -0.02 

H0.01 

-0.01 

H0.02 

'-0.1 



1.11 

1.15 



1.28: 

1.35 



In order to assess the accuracy of the IILM prediction for the nucleon at different quark masses, 
we have verified that these are consistent with the corresponding unquenched [9, 10] lattice QCD 
calculation (see Fig. 3). 
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Figure 3: The nucleon mass as a function of the pion mass squared obtained in IILM and in lattice calcu- 
lations. The blue points correspond to the IILM results, the yellow points to CP-PACS lattice result [9] and 
the red points to MILC data [10]. Since MILC and CP-PACS adopt different scheme for fixing the physical 
units, we have rescaled the lattice spacing of MILC data until the two set of lattice data are consistent. 



4. Instanton-induced dynamics in the chiral regime 

Having verified the consistency of the IILM with lattice simulations and with the behavior 
expected from #pt, we conclude that restricting to semi-classical vacuum fluctuations indeed rep- 
resents a useful approximation to the QCD path integral. Hence, we can use this model to explore 
the dynamical mechanisms controlling the role of chiral dynamics at different mass scales. An 
essential quantity in the IILM is the overlap matrix Tjj = J d 4 zYj\z)iD^YnY®(z) specifing the 
probability that a quark hops from instanton i to instanton j, this matrix plays the main role both 
in the evaluation of the fermionic determinant and of the quark propagator. In the IILM these two 
quantities can be factorized in a zero mode part and a non-zero mode part 

Det A (^ + m q ) = Det Azm x Det' Anzm , 
S(x,y) A = S zm (x,y) A + S nzm (x,y) A . 

where 

Det Azm = Det(r + m 9 ), 



S m (x,y) A = EvFtoY^OO 
u 



T — im q 



u 



The non-zero mode part of the determinant is approximated with the product of the non-zero mode 
contributions of each individual instanton: 

N + +N- 

Det Azm = n (l-34m ?P 0. (4.1) 
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Figure 4: Distribution of the maximum fermionic overlap matrix elements in the IILM for different quark 
masses. 

The non-zero part of the propagator is known analytically and is found to give sub-leading contri- 
bution to low-energy dynamics. From the structure of the quark propagator and of the fermionic 
determinant in the IILM model, it follows immediately that the role of instanton-induced chiral dy- 
namics depends crucially on the value of the typical 7y matrix element. In Fig. 4 we have plotted 
the distributions of the maximum overlap matrix elements 7)a obtained at different values of the 
quark mass. We note that such distributions are peaked around the value m* ~ 80 MeV. The value 
m* plays a central role in specifying the scale for instanton induced chiral symmetry breaking in 
the IILM. In fact, for quark masses m q S> max[7}/] ~ m* the contribution of the overlap matrix 
elements to the propagator and to the fermionic determinant becomes negligible. 

Other two interesting features of the Instanton Model can be revealed from Fig 5: (i) for m q < 
m* the density of quasi zero-modes increases as the quark mass decreases and (ii) for m q > m* the 
near zero-mode part of the spectrum becomes practically independent on the quark mass. These two 
non-trivial dynamical effects are in fact related and can be explained as follows. Quark loops are 
known to generate strong non-local correlations between pseudoparticles of opposite topological 
charge. Such correlations tend to suppress configurations in which one or more pseudoparticles 
are located far from all others. Hence, the contribution of the fermionic determinant leads to a 
reduction of the density of nearly zero-modes, as completely isolated instantons are known to have 
exact zero-modes. As the quark mass gets larger, such a topological screening becomes less and 
less effective and the population of near exact zero-modes increases, explaining the rise of the 
peak of eigenvalue density near the origin. On the other hand, we have seen that for m q > m* the 
contribution of the fermionic determinant is suppressed and fermion-induced topological screening 
completely disappears. As a result, the low-virtuality sector of the Dirac spectrum stops depending 
on the quark mass. These are a interesting predictions of the instanton liquid model that can be 
checked with unquenched lattice calculations. 

5. Conclusion 

In this talk we have used the IILM as a tool to investigate the dynamical mechanisms which 
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Figure 5: Saturation of the Dirac spectrum profile at small X for growing quark mass 



drive the dynamics and the structure of hadrons, in the light-quark sector of QCD. We have shown 
that the model's predictions are consistent with QCD in the pion mass regime < 500 MeV, where 
we expect chiral dynamics to play an important role. The model reproduces #pt prediction for the 
spectrum of the Dirac operator and is consistent with the lattice data for pion and nucleon masses 
in the range m K « 300 - 600 MeV. 

Having checked that chiral dynamics is correctly encoded in the instanton model, we have 
exploited our analytic handling on instanton-induced correlations to study the dynamical mecha- 
nisms involved in the transition into the chiral regime, at a microscopic level. We have identified 
a mass scale m* = 80 MeV above which we do not expect QCD correlators to be dominated by 
chiral dynamics and above which quenched calculations should begin to approximate full QCD. 
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